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In this work we investigate properties of fermions in the 5*0(5) theory 

of high Tc superconductivity. We show that the adiabatic time evolution 

of a S0{5) superspin vector leads to a non-Abelian SU{2) holonomy of the 

50(5) spinor states. Physically, this non-trivial holonomy arises from the 

non-zero overlap between the SDW and BCS quasi-particle states. While 

the usual Berry's phase of a S0{3) spinor is described by a Dirac magnetic 

monopole at the degeneracy point, the non-Abelian holonomy of a 5*0(5) 

spinor is described by a Yang monopole at the degeneracy point, and is deeply 

related to the existence of the second Hopf map from 5''^ to 5*^. We conclude 

this work by extending the bosonic 5*0(5) nonlinear a model to include the 

fermionic states around the gap nodes as 4 component Dirac fermions coupled 

to SU{2) gauge fields in 2-1-1 dimensions. 
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I. INTRODUCTION 



Recently, a unified theory based on S0{5) symmetry between antiferromagnetism (AF) 
and d wave superconductivity (dSC) has been proposed |jl|] for the high cuprates. Initially, 
this theory was formulated in terms of a nonlinear a model which describes the effective 
bosonic degrees of freedom below the pseudogap temperature. This theory gives a unified 
description of the high Tc phase diagram and offers a natural explanation of the vr resonance 
mode observed in the YBCO superconductors. With the exception of exact microscopic 
SO {5) models 0-0, both numerical investigations |^,^ and the experimental proposals 
Hj6| JlO|JTT[] have primarily focused on the bosonic sector of the SO (5) theory. 

However, it is clear that a complete theory of high Tc superconductivity has to properly 
account for the fermionic degrees of freedom as well. Some key experiments on the pseudogap 
physics, e.g. the ARPES experiments, primarily probe the single electron properties rather 
than the collective modes. Within the 5*0(5) theory, the pseudogap regime is identified with 
the fluctuations of the orientation of the 5*0 (5) superspin vector. Therefore, it is essential to 
understand how the fluctuations of the superspin couple to single particle fermionic degrees 
of freedom. Because of the d wave nodes, there are fermionic excitations with low energy, and 
they make important contributions to thermodynamics and to the damping of the collective 
modes. 

Motivated by these considerations, we investigate the fermionic sector of the 5*0(5) 
theory. Our main interest is to understand how the fermionic degrees of freedom are coupled 
to the bosonic 5*0(5) superspin vector, and investigate the novel topological properties 
of this coupling. Rather surprisingly, we find that this coupling leads to a SU{2) non- 
Abelian holonomy of the fermionic states. The fermionic states in the 50(5) theory are 
nothing but the familiar SDW (spin-density-wave) and BCS (Bardeen-Cooper-Schrieffer) 
quasiparticles relevant for the AF insulator phase and the dSC phase. Since AF and dSC 
states have very different physical properties, one would naively expect the SDW and BCS 
quasiparticle states to be orthogonal to each other. The central result of our work shows 
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that this is not the case. The SDW and BCS quasi-particle states have nonzero overlaps, and 
this overlap defines a precise connection for their adiabatic evolution. This simple physical 
property of non-orthogonality of the SDW and BCS quasi-particles leads to an extremely 
rich mathematical structure. In particular, because of the inherent degeneracy of the SDW 
and BCS quasi-particle states, the adiabatic evolution is non-Abelian, and can interchange 
the degenerate states upon a cyclic evolution returning to the origin. 

Our current investigation is also motivated by the question "what is special about the 
SO {5) symmetry group?". Is it simply introduced as a convenient mathematical description 
of the AF and dSC phases in a unified framework, or is there something deeper which explains 
its uniqueness and calls for its natural emergence? In this work, we give partial mathematical 
answers to these probing questions, as we embark on a journey through some of the most 
elegant and beautiful mathematical concepts in group theory, differential geometry, topology 
and division algebra, which, as shall see, are unified by the SO (5) group in a profound and 
unique way. The following table summarizes some of the main mathematical properties of 
the holonomy of a SO {5) spinor, in comparison with the Berry's phase of a 5*0(3) spinor: 





SO (3) Spinor 


50(5) Spinor 


Holonomy 


U{1) Berry's Phase 


SU(2) Wilczek-Zee Holonomy 


Singularity at Degeneracy Point 


Dirac Monopole 


Yang Monopole 


Topological Invariant 


First Chern Number 


Second Chern Number 


Wigner-von Neumann Class 


Unitary 


Symplectic 


Associated Hopf Maps 


53 to 52 


S^ to S^ 


Associated Division Algebra 


Complex Numbers 


Hamilton Numbers (Quaternions) 



Soon after Berry's discovery of the adiabatic phase [jT2|,|T3[, Wilczek and Zee general- 
ized this concept to the non-Abelian holonomy in a quantum system with degeneracy. In a 
sense, the most natural generalization of the concept of the Berry's phase of a 5*0(3) spinor 
is the SU{2) non-Abelian holonomy of a 5*0(5) spinor. The usual Abelian Berry's phase 
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has its mathematical origin in the first Hopf map from a three sphere 5*^ to a two sphere S"^. 
Similarly, the non-trivial holonomy of a 5*0(5) spinor is deeply related to the existence of the 
second Hopf map from the seven sphere S*^ to the four sphere S*^, the later being the order 
parameter space of the 5*0(5) theory. The generalization from the first Hopf map —>■ S*^ 
to the second Hopf map S*^ S*^ is uniquely related to the generalization from complex 
numbers to Hamilton numbers (quaternions). In mathematics, there exist only three kinds 
of division algebra, namely complex numbers, Hamilton numbers (quaternions) and Cayley 
numbers (octernions). These three classes of division algebra lead to three types of Hopf 
maps |jl5|, —>■ 5*^, 5*'' —>■ S*^ and S^^ —>■ S*^. Octernions and the third Hopf map are not 
useful in physics because of the lack of associativity. This makes the second Hopf map and 
its associated non-Abelian holonomy of a 5*0(5) spinor an essentially unique generalization 
of the usual Abelian Berry's phase. 

The novel topological property of a 5*0(5) spinor enables us to extend the bosonic 5*0(5) 
non-linear a model to include the fermionic states. Because of the d wave nodes, the low 
energy electronic states can be best described as Dirac fermions in 2 + 1 dimension. The 
SU (2) non-Abelian holonomy uniquely determines the coupling of these Dirac fermions to 
the fluctuation of the 5*0(5) superspin order parameter, and takes the form of a minimal 
coupling to a SU{2) gauge field. With the inclusion of low energy fermionic modes, the 
formulation of the 50(5) theory is essentially complete. The present formalism can be 
used to systematically investigate the fermionic properties which result from the interplay 
between AF and dSC, for example, the pseudogap physics, fermionic quantum numbers and 
bound state inside the 5*0(5) vortices and junctions, electronic states in the stripe phase 
etc. The non-Abelian holonomy of the BCS and SDW quasi-particles discussed in this work 
could also be observed directly in experiments. However, in the paper, we shall only restrict 
ourselves to the mathematical formulation of the theory, physical application of the present 
formalism will be discussed in future works. 

Some of the mathematical properties presented in this paper have been discussed previ- 
ously in other contexts. General background on Berry's phase and its various generalizations 
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have been collected in an authorative reprint volume edited by Shapere and Wilczek [jT6| 



Minami studied the connection between the second Hopf map and the Yang monopole. 
Wu and Zee [ff8[ studied a 5*0(5) nonlinear a model in seven dimensions and discussed the 



connection between the second Hopf map and quaternionic algebra. Avron et al |T9|J20 
discussed the general settings of Berry's phase in a fermionic time reversal invariant sys- 
tems and the connections to Wigner von Neumann classes. An interesting example of the 
the non-Abelian Berry's phase has been studied by Mathur, who showed that such term 
appears in the adiabatic effective Hamiltonian for the orbital motion of a Dirac electron 
and leads to spin-orbit interaction and Thomas precession ||2l|]. Shankar and Mathur later 



identified the non-Abelian Berry vector potential in this problem with that of a meron [22 



While some of the mathematical concepts discussed in this paper may not sound familiar 
to some readers, we shall present our paper in an essentially self-contained fashion which 
should be accessible without much background mathematical knowledge. 



II. HOLONOMY OF A SO(5) SPINOR 

A. Adiabatic Evolution of a 50(5) Spinor 

Abelian Berry's phase or f/(l) holonomy is a familiar object in physics. A typical ex- 
ample comes from considering a spin ^ particle coupled to a magnetic field p!^JT3|j2^ . The 
Hamiltonian for this problem is a simple 2x2 matrix, given by 

Ha/3 = A na{t)alp (1) 

where na{t) (a = x,y,z) is a three dimensional unit vector, a" are the three 2x2 Pauli 
matrices and A is the Zeeman energy gap. The projection operator 

P± = l{l±n^cr-) (2) 

projects onto the subspace of non- degenerate eigenvalues ±A. Under the assumptions of adi- 
abatic evolution, we can choose the instantaneous eigenstate corresponding to the eigenvalue 
+A as 
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(3) 



where Ni is a normalization factor, which ensures (^V^^(t)I^V^^(t)) = 1. The Berry's phase 
is defined in terms of the non-vanishing overlap between these instantaneous eigenstates: 

d 



l?{t) = -^l rfr(vl/«(r)|-|v[/«(r)) 



(4) 



This phase factor can in turn be expressed as a line integral over the vector potential of a 
Dirac monopole: 



7i'^(t) = / dTAl{n{T))dnaiT) 
Jo 



(5) 



where the vector potential 



Aldua 



1 



-{n^duy - Uydn^) 



(6) 



2(1 + n,) 

is singular at the south pole = —1. 

Similarly, one can choose another class of instantaneous eigenstates corresponding to the 
+A eigenvalue, defined as 



1 



(7) 



The Berry's phase associated with this class of instantaneous eigenstates is given by 
(2) / . /■* . , ^ d 
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\t) = -il dT{^'^l\T)\ — \^f{T))= I dTAl{n{T))dnJ 
Jo dr Jo 



T 



(8) 



where the vector potential 

Aldua 



1 



{n^dny - Uydn^) 



(9) 



2(1 

is singular at the north pole ^ = 0. A^{n) and A^{n) define the two nonsingular "patches" 
of the monopole section in the sense of Wu and Yang ||25|, their difference is a pure gauge 
in the overlapping equatorial region S^, where = and n"^ + riy = 1: 
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dW = A'^dua — A\dna = Uxduy — nydn^ (10) 

Obviously, we can also define the Berry's phase of the instantaneous eigenstates with 
eigenvalue —A, and find their associated gauge potentials A_{n) and A_{n). They corre- 
spond to a Dirac monopole with opposite magnetic charge. 

How does the concept of a Berry's phase of a 5*0(3) spinor generalize to the case of 
a SO {5) spinor? The Hamiltonian for a S0{5) spinor coupled to a 5*0(5) unit vector 
(a = 1, .., 5), called superspin in the 50(5) theory, is given by a simple generalization of (|lj): 

n^p = A n,it)Tlp (11) 

where the five 4x4 Dirac F matrices are given by 







-lOy \ 


p(2,3,4) ^ 


a 

























(12) 



Here a = (cr^, ay^a^) are the usual Pauli matrices and *(? denotes their transposition. A has 
the physical interpretation of a SDW gap energy when the superspin vector points in the 
2, 3, 4 direction, and a BCS gap energy when it points in the 1, 5 direction. 

Here we see a crucial difference between the 5*0(3) spinor Hamiltonian (|I]) and the 5*0(5) 
spinor Hamiltonian (|ll|). The instantaneous eigenvalues of both Hamiltonians are ±A. 
However, the eigenvalues are non-degenerate for the 50(3) case, but doubly degenerate for 
the 5*0(5) case. For example, when the n-field points along n^, Hamiltonian (pJ]) is diagonal, 
with Till = 7-^33 = —7^22 = —^^44 = A. It is then immediately obvious that states in the 
subspace of |\l/2) and |\E'4) have energy —A, and states in the subspace of and |\l/3) have 
energy +A. 

Wilczek and Zee |T^ pointed out that adiabatic Hamiltonians with degeneracies can have 
non-Abelian holonomies and may be characterized by the non-Abelian gauge connections. 
Let us imagine a slow (adiabatic) rotation of the 50(5) unit vector field when it traces a 
closed cycle on 5*^. What happens to our our spinor states? In the adiabatic approximation 
we can assume that when we start with a spinor state in the low energy subspace it will never 

7 



be excited into the upper subspace. However the low energy subspace is two-dimensional. In 
the course of adiabatic rotation, the spinor state is always some linear combination of the two 
degenerate low energy states. Different instantaneous eigenstates are related to each other 
by a 2 X 2 unitary matrix. Moreover, after the five-vector returns to its original direction 
the spinor state does not necessarily return into itself but into some linear combination of 
the initial state with its degenerate orthogonal state. Therefore spinors in SO (5) have a 
non-Abehan SU (2) holonomy which has its origin in the double degeneracy of spinor states 
and leads to many novel phenomena that we discuss below. 

Following the discussions of the Abelian Berry's phase, we introduce the projection 
operators 

p± = i(i±n„n ; Pl^P± ; np±^±AP± (13) 

It is convenient to parameterize the four sphere as 

rii — sin9isin92sin(f)2 
n2 — sin9icos92Cos(f)i 

— sin9icos92sin(pi 
rii — cos9i 

— sin9isin92Cos(f)2 (14) 

With these representations, we can choose one class of instantaneous two dimensional basis 
states corresponding to the -|-A eigenvalue as 

1 T 

^=^= (l + cos9i, sin9icos92e''''\0, isin9isin92e'^A 
2(1 + cos9i) ^ ' 



1 T 

(O, isin9xsin92e-'"*"' , 1 + cos9x, sin9isin92e~'^') (15) 



^2(1 + cos9i) 
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where we defined unit vectors such that {ei)j — 5ij (for example ei = (1,0,0,0)-^) and 
normahzation factors A^i and A^3 are chosen such that (^^''l^^-') = 6ij. These states obey 
H{n)\'^^^^) = +A 1^'+'') by construction. From this basis of instantaneous eigenstates we 
obtain the following non-Abelian holonomy matrix: 

Alternatively, we could have used the following set of instantaneous eigenstate basis for 
the +A eigenvalue, 



iA\dna 



(16) 



Prom this alternative basis set we obtain the following non-Abelian holonomy matrix: 



(17) 



(18) 



B. Adiabatic Connection and Yang's Monopole 



Prom the discussions in the previous section, we see that the non-Abelian vector potential 
A"' and A"' are direct generalizations of the vector potentials ^4" and in the Abehan Berry's 
phase problem, which are the vector potentials of a Dirac magnetic monopole. Furthermore, 
the 2x2 matrix A"'dna is traceless and hermitian, thus defining a SU (2) gauge potential. 
Therefore, it is natural to ask if A"' and A"' are the SU{2) generalizations of the Abehan 
Dirac monopole potential. 

Direct computation of A"' gives 



i[Al]udna = 



2(1-1- sinOcosa) 



i[Al]i2dna 



2(1 + 714) 

^g-*(</'l+<A2) 

2(1 + 0056*1) 



{n2dnz — n^dn2 — riidn^ + n^drii) 



{sin^Oi d92 + isin^9isin92Cos92 d(pi — isin^9isin92Cos92 dcpo^ 



-{n2dn^ — n^dn2 + riidn^ — n^drii 



2(1 + ^4) 

+i {nidn2 — n2dni + n^dus — n^dn^)) 

i[Al]2idna = i[Al]l2dna 

t[Al]22dna = -i[Al]ndna (19) 

The explicit form of A^ can also be obtained and it may be shown that it has singularity 
at ?T,4 = 1. Therefore, while A+ is non-singular except at the "south pole" ^4 = —1, A+ is 
non-singular except at the "north pole" n4 = 1. In the overlapping region, they are related 
to each other by a pure non-Abelian gauge transformation: 



A^dua = W-\Aldna) W - t W-^dW 

1 - ins Til - in^ 

1 2 

1 — 77.4 —m — n2 + in^ 



(20) 



Similarly, we can define the vector potentials A-. and A- associated with the —A eigenvalue. 
These gauge fields turn out to be exactly the vector potential of a Yang monopole! 

In 1978, Yang found a beautiful generalization of the concept of the Dirac magnetic 
monopole p6| , p7| . While the Dirac monopole is a S0{3) symmetric point singularity in the 



three dimensional space, the Yang monopole is a 5*0(5) symmetric point singularity in the 
five dimensional space. The Dirac monopole defines a topologically nontrivial U{1) fiber 
bundle over the two sphere 5*^, the Yang monopole defines a topologically nontrivial SU{2) 
fiber bundle over the four sphere S*^. The best way to understand a magnetic monopole 
is through the concept of Wu-Yang section. For the case of a Dirac monopole, one divides 
the S"^ sphere into a "northern hemisphere" and a "southern hemisphere", over which two 
different non-singular vector potentials A"' and A"' are defined. The "overlapping region" 
between the two hemispheres is the equator, or 5*^. In this overlapping region, the two 
f/(l) gauge potentials can be "patched together" non-trivially A"- = A"- — i W'^&^W . Since 
the mapping from the overlapping region S*^ to the group space U{1) can in general be 
characterized by a winding number, this integer uniquely defines the non-trivial U{1) fiber 
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bundle over S'^. This integer is called the first Chern number, and is defined by 



27rci 



i 
i 



a 



ab 



i 



A^dn, 



'a 



= I 



i 



W-^d^'Wdn, 



a 



(21) 



where F'^^ = d'^A^ — d^A'^ = d°'A^ — d^A"- is the field strength, W is a gauge transfor- 
mation between A and A defined in equation (0) and daab is a surface element of a 2- 
sphere. Here and everywhere else in this paper we normalize surface element on a n-sphere 
S"" by requiring that when integrated over a unit sphere it gives the correct surface area 
is„ e'^i -""+iXai(i(Ta2...a„+i = = 27r("+^)/^/r( ^^) . Rcadcrs familiar with differential geome- 
try will easily recognize daa2...an+i = dxa2 A . . . Adxa„+i- For the Berry's connection of a spin 
one half particle, we obtain from formula (^) ci = ±1 for the A"^ and A"i vector potentials 
respectively. 

From the point of view of Wu-Yang section, it is straightforward to see why a SU (2) 
point monopole has to be enclosed by 5*^. If we view the "north pole" as 724 = +1 and 
the "south pole" as 71,4 = —1, then the "overlapping region" between the two hemispheres is 
given by ?T,4 = 0, and is a three sphere defined by n^ + n^ + ^l + ng = 1. Since the mapping 
from to the group manifold of SU (2) can be characterized by a winding number, one 
can define a topologically non-trivial "patching" between the "northern hemisphere" SU(2) 
gauge potential A and the "southern hemisphere" SU (2) gauge potential A. This integer is 
called the second Chern number and is defined by the generalization of (pl|): 



where J-'"''^ = d"'A^ — d^'A'^ + is the non-Abelian field strength associate with the 

vector potential A" or A" [g2|,0. 




(22) 
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From the above equation, we easily recognize the integral over the three sphere as the 
non-Abelian Chern-Simons term. One may wonder in what sense is the Chern-Simons term 
a natural generalization of the Bohm-Aharonov type of line integral defined in (pT]). To make 
the connection between (|21| ) and ( p2D more precise, one can introduce an Abelian totally 
antisymmetric three-index gauge field defined by 

3 

= 1 ^ (_)^ Tr(^^('^)5^(^)^^(^) - - (23) 

2 P{abc) 2 

and its associated four-index field strength F"'^^'^. In this representation, the second Chern 
number appears to be a direct extension of the first Chern number: 



87r^C2 = i A'^^'daabc - f A'^'^da,,, = i F'^''^" da abed (24) 
Js-i Js^ Js-^ 

Since the SU (2) holonomy connection we found explicitly in ([T9| ) is 5*0(5) symmetric and 

gives C2 = ±1 for the A+ and A^ gauge potentials respectively, it can be uniquely identified 

with the gauge potential of a Yang monopole. The introduction of the Yang monopole into 

the 50(5) theory is a important step in describing the fermionic degrees of freedom. A 

single bosonic 5*0(5) rotor has only the fully symmetric traceless tensor representations of 

the 5*0(5) group. However, a single rotor with a Yang monopole at the center contain all 

irreducible representations of the 50(5) group, including the fermionic spinor sector |27||. 



C. Spinor Rotation Matrix 

Having established the topological structure of the non-Abelian holonomy, let us now go 
back to our problem of the coupled 5*0(5) spinor and vector and try to actually solve the 
time-dependent Hamiltonian 

Hit) =Ana miKf,^^ (25) 

For the purpose of later application to BCS and SDW quasi-particles, we use here the 
second-quantized notations, so \E'q,'s are now operators that may be thought of as creating 
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spinor states out of the vacuum. Time dependence of these operators may be found from 
the Heisenberg equation of motion by taking a commutator of \I/q, with the Hamiltonian and 
we shall assume that \l/a's anticommute with each other. We decompose \E'a using 

^a(t) = Sa(sm(s{t) (26) 

where matrix Saf3{t) is a unitary matrix and is another spinor. The Heisenberg equation 
of motion for \1/q, can be expressed as 

Sa(3^f3 + Sa0^(3 = A UaTl^Sp^^^ (27) 

If we choose such that its time evolution is given by the time independent Hamiltonian 

no = A $t r^,$^ (28) 

then ( P7| ) can be solved if Sa/s satisfies the following two conditions 

n,(t)St(t)r'^S(t) = (29) 
S^{t)dtS{t) = (30) 

The meaning of decomposition ( PB| ) is clear. At each moment t we use matrix Sap to 
rotate the spinor from its instantaneous direction na{t) to a fixed direction n4. However this 
does not define the matrix Sap uniquely. There are two states pointing along 714 and two 
orthogonal states pointing in the opposite direction. So the first equation ( ^9]) only fixes Sap 
up to arbitrary rotation within these two pairs. The second equation (^) combined with 
the first one determines Sapit) uniquely. 

Let us now see how we can construct the spinor rotation matrix that satisfies both 
conditions (^), (pOl). For a moment we forget about the second condition and try to 
find any unitary matrix Sap such that ( pQ]) is satisfied. Matrix has two eigenvectors 
with eigenvalues +1 and two eigenvectors with eigenvalues -1. Therefore if we find two 
orthogonal eigenvectors of n^F" with eigenvalues +1 and two with eigenvalues -1, they will 
specify a necessary rotation for us ( eigenvectors that correspond to different eigenvalues are 
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orthogonal ). These eigenvectors may be easily found using projection operators P± defined 
indl). 

We take P+ei , P+e^ , P„e2 , P^e^ and call the normalized vectors ?/'•*•'. They satisfy 



yii)*yU) = A. 



where Aj = +1 when i = 1,3 and Aj 



(31) 



-1 when i = 2,4. It is then clear that we should 



define Sap = y^^- For example, if we parameterize superspin as in ([l4|) , the matrix Sa/3 is 
explicitly given by 

—sin^-cos02e~^'^^ isin^-sin02e~'^'^'- 



a/3 



2 — 2<: 

isin^-sin92e^'^^ 



cos- 



sin^cos92e~ 



COS 



01 



(32) 



Now we recall about the second condition (|30| ) on 5*0,^. Here it is important to keep in 
mind that we are working within adiabatic approximation where transitions between states 
of different energies are forbidden. Therefore in equation ( pO|) we don't need to consider the 
full 4x4 equation but only two 2x2 equations (one of the positive energy subspace with 
components 1 and 3 and the other of the negative energy subspace with components 2 and 
4). For notational convenience let us introduce a matrix 



G 



V 



1 



(33) 



^ 1 ^ 
10 
10 
1 

From Section |111 A| we know that within positive and negative energy subspaces S'^dS cor- 
responds to an infinitesimal SU{2) rotation, therefore, it can be expressed as 

^ idA+ ^ 
idA- 



SUS = G 



V 



G 



(34) 
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where the 2x2 SU{2) gauge matrices dAj^ and dA- are specified in section [III A| . Using S", 
w4+, and A- we can finally express S as 



Sit) = S{t)G 



1+ 
W_ 



a (35) 



where 

dU±{t) = -% U±{t)dA± ; U±{t) = Te-'/o ■^iW'^"'"^ (36) 

T is the time ordering symbol and matrices U± are the finite SU{2) holonomy matrices for 
the upper energy and lower energy subspaces. It is easy to check that the resulting S{t) 
satisfies both equations (^) and (|30|) . 



D. Hamilton's Number, Hopf's Map 

In 1843, while walking along the Bridge of Brougham, Hamilton discovered a beautiful 
generalization from complex numbers to quaternions; In 1931, the same year when Dirac 
formulated the theory of the magnetic monopoles, Hopf introduced the concept of Hopf 
maps. These two great discoveries in algebra and geometry are intimately related to each 
other [|l^,|l^,|l^ , and they provide a deep explanation for the non-Abelian holonomy of a 



5*0(5) spinor. 

The first Hopf map is a mapping from 5''^ to S"^ and is related to Dirac's magnetic 
monopole, while the second Hopf map is a mapping from 5''' to S^, and is deeply related 
to Yang's monopole. is locally isomorphic to S*^ x f/(l), where we can view S"^ as the 
sphere enclosing the Dirac monopole and U{1) is the gauge field due to a Dirac monopole. 
In the presence of a Dirac monopole, the U{1) bundle over S*^ is topologically non-trivial. 
However since is parallisible, one can use the first Hopf map to define a nonsingular vector 
potential due to Dirac's monopole everywhere on S^. Similarly, 5''' is locally isomorphic to 

X SU{2), and one can use the second Hopf map to define a non-singular SU{2) gauge 
field everywhere on S'''. 
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It is simple to define the first Hopf map in physicist's language. One introduces a two 
component complex scalar field with a constraint, 



Then the first Hopf map is defined as 



zh=\z^\^ + \z2Y = l (37) 



?7i : na = z^OaZ (38) 

where ex's are the usual Pauli matrices. Condition (|37D ensures that r? = 1, so we have a 
mapping from z with three independent components ( ) into n"- with two independent 
components { S"^ ) rji : ^ S^. 

If we start from the singular gauge potential on S'^ defined by (H) and substitute for 
Ua the definition of the Hopf map (P3), we find an induced gauge potential on 



1 ^2 

ui = iA'^dua = uj^dza + ujadz]^ = —-\z2dz2 — Z2dz2 — -. — 7z{z\dzi — Zidz\)] (39) 

2 I Zip 

This U{1) gauge potential is singular on a great circle \zi\ = of S^. However, unlike 
its counterpart A" on 5*^ , the singularity of uj can be completely removed by a gauge 
transformation: 

Zl 1 
M = 1 — r ; n = uj + u^^du = —-{zj^dza — Zadzl) (40) 

We see that the new U{1) gauge potential ft is non-singular everywhere on S^. Similarly, 
the gauge potential lj induced by A" defined in (|]) is singular on the great circle \z2\ = of 
S^. This singularity can also be removed by a gauge transformation using u = f^, giving 

I ^2 I 

the same non-singular U{1) gauge potential fl. 

This simple calculation demonstrates the deep relation between the first Hopf map, 
Dirac's monopole and the Berry's phase. We can think of the first Hopf map as defining a 
relation between a S0{3) vector Ua and a 5*0(3) spinor {zi,Z2). This relation is invariant 
under a U{1) gauge transformation Za — > e^'^Za- If the Ua vector goes from point A to point 
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B on S'^ via two different paths, the Za vector is transported from point a to two different 
points bi and 62 on S^. (See figure |1|). 




FIG. 1. Transport of a vector on S'^ defines transport of a spinor on S^. The final spinor state, 
i.e. the final point on depends not only the final point of on but on the path itself. Points 
61 and 62 on have the same image B on 5^ under the first Hopf map. 

The transport of Za is uniquely defined by the non-singular gauge potential Q. Since 
the two points 61 and 62 differ from each other by a pure phase, they are projected to the 
same point B on S'^. Their phase difference is exactly the Berry's phase. Since 5*^ is locally 
isomorphic to S"^ x U{1), it encapsulates the full quantum information, including both the 
direction of the SO (3) vector and the phase of the 5*0 (3) spinor. This way, we can also 
understand the singularity associated with a Dirac monopole as arising from the projection 
from a non-singular U{1) gauge field on S^. 

There is an elegant generalization from these considerations to the case of a 5*0(5) spinor. 
The SU{2) holonomy of a SO (5) spinor and the Yang monopole are deeply related to the 
second Hopf map from S*^ to S"^ and a non-commutative division algebra called quaternions. 
A quaternion with three imaginary units i, j and k is a generalization of a complex number 
with one imaginary unit i. These three imaginary units mutually anticommute and they 
obey the following multiplication table: 

i2=f = k2 = -l ; ij = k ; jk = i ; ki = j (41) 

It is simple to understand quaternions from a physicist's view point, since the algebra of 
these three imaginary units S'^ = (i,j,k) is identical to the algebra of three Pauli matrices 
cr'*', and a one-to-one correspondence between them is therefore possible. A quaternion can 
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be expressed as 



q = + i + j + k 



(42) 



and therefore has four real components. Complex conjugation of quaternions is defined as 
the operation that changes the sign of the imaginary parts, i.e. of the last three components. 
Magnitude of the quaternion is given by |qp = q*q. Using a two component quaternion, we 
can parameterize the S"^ sphere by 



qo 



q^q = |qiP + |q2p = i 



(43) 



With this notation, the second Hopf map can be expressed as a simple generalization of 



112 ■ na = qSaq 



(44) 



Here is a five dimensional real vector, and the five 2x2 quaternionic valued 7 matrices 
are defined as a simple generalization of the Pauli matrices 
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From the norm of qo, defined in 
show that = 1. Therefore 



(45) 



I) and the properties of the 7 matrices, one can easily 
defines a mapping from S"^ to S^. It may be compared 



with Section |11 A| where an S0{5) spinor was represented by a four component complex 



vector and equation (0) implied that each spinor was mapped into an 5*0(5) vector by 

Ha = ^^r,^ (46) 

Connection between two representations is easily established by choosing 

^ + ,^2^ ^1 + ,^2^ qO _ ^^3^ ^0 _ ^g3| (47) 



which makes (^4]) and (^6]) identical, and turns the second equation in (^31) into a normal- 
ization condition for the spinor wavefunction. 

Any SU{2) gauge field A can be expressed in terms of the Pauli matrices as A"" = 
A"'^a\. Because of the isomorphism between the Pauli matrices and the three imaginary 
quaternionic units S'^ = (i,j,k), a SU{2) gauge field can also be expressed as a imaginary 
quaternionic field: 

i A" = A^^Sx (48) 

Using this observation and directly substituting the definition of the second Hopf map (^i]) 
into the SU (2) gauge field of a Yang monopole (pJf), one finds that it induces a SU (2) gauge 
field on S*^ defined by 

uj = iA^dua = ujldqa + ujadql = - |2 [qi(q2^q2 - '^q2q2)qi - |q2p(c?qiqi - qic^ql)] (49) 

^|qi I 

where we established correspondence between quaternions and Pauli matrices as 

i icTy j -ian, k ^ -ia^ (50) 

The SU{2) gauge field in (^) is singular on a three dimensional "great sphere" |qi| = 
of S*^. Fortunately, just like the case of the first Hopf map, the singularity can be completely 
removed by a SU{2) gauge transformation defined by a unitary quaternion u = qi/|qi|. 
The resulting SU{2) gauge field 

Q, = u'^LJU + u^du = {qidqa - «^qLq«) (51) 

This SU{2) gauge field over S'^ appears as a direct quaternionic generalization of the U{1) 
gauge field over S"^, and it is non-singular everywhere. Similarly, one can map the singular 
SU{2) gauge field A+ on to a singular SU{2) gauge field Cj on 5^. Upon removing the 
singularity by a unitary quaternion u = q2/|q2|, the resulting SU(2) gauge field is again 
given by the non-singular gauge field f2. 

Like the previous calculation, this calculation demonstrates the deep relation between 
the second Hopf map, Yang's monopole and the SU{2) holonomy of a S0{5) spinor. We 
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can think of the second Hopf map as defining a relation between a 5*0(5) vector and 
a normahzed 50(5) spinor = (qi,q2)- This relation is invariant under a quaternionic 
unitary transformation qc q^u with u"''u = 1. Such a quaternionic unitary transformation 
is identical to a SU (2) gauge transformation. In this case, if the Ua vector goes from point 
A to point B on via two different paths, the q^ vector is transported from point a to two 
different points hi and 62 on S*^. (See figure |]). 

Yang Monopole 



s 

FIG. 2. When a vector is transported on along two different paths a spinor is transported on 
5*^ to two different final states bi and 62, which have the same image B on 5'^. 

The transport of qa is uniquely defined by the non-singular gauge potential fi. Since 
the two points hi and 62 differ from each other by a SU (2) gauge transformation, they are 
projected to the same point B on S*^. Their difference is exactly the SU{2) holonomy of a 
50(5) spinor. Since 5''' is locally isomorphic to S"^ x SU (2), it encapsulates the full quantum 
information, including both the direction of the 5*0(5) vector and the SU{2) holonomy of a 
5*0(5) spinor. 

Another way of seeing the invariance of the second Hopf map under SU(2) gauge trans- 
formation comes from returning to a representation of the 50(5) spinor as a four component 
complex vector ^. Choosing components of as 

^ = {iJi,(Pi,4,4f (52) 
gives for the second Hopf map Ua = "^^Ta"^ 

ni = (j)itjJa + ^l0Q 

n2 = i{4>iiJa - i^i<Pa) 
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^5 = (Pi(pa - ^i^a (53) 

One can easily see that this exhausts all possible bilinear combinations that are both real 
and invariant under simultaneous SU(2) rotations of and ip [ p8|] . 

To conclude this section we remark that analysis presented here leads us to an intriguing 
conclusion that not all spheres are created equal. and S"^ are endowed with a unique 
algebraic structure and associated non-singular complex (140|) and quaternionic (|5l|) gauge 



connections. This deep connection between the algebra and geometry underlies the non- 
trivial holonomy of SO {3) and 5*0(5) spinors. 

E. HP{1) Representation of the SO(5) Nonlinear a Model 

It has long been known that the 5*0(3) nonlinear sigma model used to describe the low 
energy physics of magnetic systems has an interesting representation due to the existence 
of the first Hopf map: OP(l) representation p9|. Indeed one observes that when we define 



the first Hopf map as in equations ( pTj) and ( |38D of Section |11 D| , the partition function for 



the nonlinear sigma model may be written as a path integral over z and gauge field A 

This may be easily proven by noticing that the action on the right hand side of ( ^41) is 
quadratic in A, so integration over the gauge field amounts to taking the saddle point value 



M 2 



z' 



^d^z - {d^z^)z] (55) 



After that simple algebra gives 

1 



^i^,nr = \i^,-^A,)z\' (56) 
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which proves equation (^41). Notice that the gauge field in (|55D is the same as in which 
is not surprising since the origin of both gauge fields is the invariance of (^Sj) and ( |56| ) under 
gauge transformations 



A^-^A^ + d^(j) (57) 
Analogously one can prove that the second Hopf map, defined by p3|) and in Section 



UP , gives rise to an HP{1) representation of the 5*0(5) nonlinear sigma model 

D^q = d^ci-qB^ (58) 
where similarly to CP{1), the HP{1) action 

i(9X)' = iD,q)HD,^) (59) 
as well as the second Hopf map (|^) are invariant under gauge transformation 

^ u-^B^u + u-'d^u (60) 

with u being any quaternion that satisfies u^u = 1. There is, however, an important 
difference between gauge transformations in (^) and (^Tf). In the former case they are 
non-Abelian SU(2) gauge transformations whereas in the latter case they are Abelian U{1) 
gauge transformations. 

It is also not surprising that the saddle point value of the gauge field B in (0) is given 

by 

B/. = [ (lid^fla - {d^qL) ] (61) 



which agrees exactly with expression ( |5TD that we obtained in the study of the 5*0(5) spinor 
holonomy (notice that both B and the one form O in ( pT| ) are quaternionic imaginary, or 
anti-hermitian) . 
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So far the HP{1) model is only an interesting reformulation of the 5*0 (5) nonlinear a 
model. However, it has the advantage that possible topological terms can be easily identified 
within this formulation. In 2 + 1 dimension, one could add the second Hopf invariant, or 
the non-Abelian Chern-Simons term 

£^TrJ (fx e^^" (b^9,Ba - ^B^B.Ba) (62) 

to the HP{1) action (^8l), where Tr implies taking the real part of a quaternion ( in repre- 
sentation where quaternions are 2x2 matrices as in (|50|) , it corresponds to actually taking 
a trace ). Because the Chern-Simons term changes by 

^TrJ d'x t^^^ {u-\d,u)u-\d,v.)vr\dxu)) (63) 

under global gauge transformation (0), the coupling constant n has to be quantized ( for 
n = 1 expression in ( |63D is a winding number of the gauge transformation which is always 
integer), i.e. only integer values of n are allowed. In 3 + 1 dimension, one could add a 9 
term, 

d Tr J d'x e^^^'^F^.FA. 
F^, = 9^B, - d,B^ + [B^, B,] (64) 



to the action HP{1) action (pSf). In 6 + 1 dimensions, Wu and Zee [|T^] showed that one 
could also add an Abelian Chern-Simons term for the 3 index gauge field introduced in 
(p3|). However, we haven't yet identified microscopic models which would give rise to these 
topological terms in the action. 



F. Wigner - von Neumann Class 

The non-trivial holonomy of the quantum mechanical wave functions are intimately re- 
lated to the degeneracy point of the Hamiltonian. Wigner and von Neumann classified the 
degeneracy points into three classes according to the generic symmetries of the Hamilto- 
nian. Time-reversal invariant systems without Kramers degeneracy (system of bosons or 
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even number of fermions) belong to the orthogonal class. The generic degeneracy point has 
co-dimension 2, i.e. one needs to tune two parameters simultaneously to zero in order to 
reach a point of degeneracy. Time-reversal breaking system belong to the unitary class, 
and the degeneracy point has codimension 3. Berry's phase of a 5*0(3) spinor describes the 
non-trivial holonomy around this point singularity in the 3 dimensional parameter space, 
which can be described as a Dirac monopole. Time-reversal invariant system with Kramers 
degeneracy (e.g. system with odd number of fermions) belong to symplectic class. The 
degeneracy point for this class has codimension 5. Unlike the two previous classes, we have 
to consider a 4 x 4 matrix problem in this case, in order to describe the level crossings of 
two different Kramers doublets. In his seminal work on random matrix theory, Dyson ||30 



noticed that this situation can be best described as a 2 x 2 quaternionic matrix problem. In a 



thorough analysis, Avron et al |[T9|j20|l showed that the natural symmetry for the symplectic 
class is therefore a unitary rotation of a two component quaternion, which is isomorphic 
to a 50(5) symmetry group. The non-trivial holonomy of a 50(5) spinor arises when one 
encircles this degeneracy point in 5 dimensional space, which acts as a Yang monopole. 

In this context, it is highly remarkable that 5*0(5) symmetry emerges naturally in a 
generic fermionic system with time reversal symmetry. To make the above discussion more 
explicit, let us recall that close to the degeneracy point, a generic Hamiltonian without 



Kramers degeneracy can be expressed as a 2 x 2 Hermitian matrix 

a, ben, d* =ceC (65) 



^ a d^ 



n 



whose eigenvalues are 



a + b±^{a + b)'-4{ab-m 
A± = (66) 

Degeneracy requires (a — 6)^ + 4|(ip = 0, which is only satisfied when a = b and Re d = 
Im (i = 0. When the Hamiltonian is real ( for systems without time- reversal symmetry 
breaking ) we have only two conditions for the degeneracy point, i.e. the co-dimension 
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is 2. But when the time reversal-breaking is present, we have three conditions and the 
CO dimension is 3. 

Up to a overall additive constant, this 2x2 Hamiltonian can be mapped onto a effective 
spin 1/2 problem in a magnetic field: Tiafi = na<Jai3, where Ua labels the three directions 
away from the degeneracy point = 0. Wave functions defined for parameters on S"^ 
enclosing the singularity are non-degenerate, but experience non-trivial holonomy due to 



the degeneracy point, as discussed in Section [II A . 



Close to the degeneracy point of two Kramers doublets, the Hamiltonian is a 4 x 4 matrix, 
and can be represented by the matrix of a spin 3/2 particle. However, because of the time 
reversal symmetry, the Hamiltonian has to be quadratic in these spin 3/2 matrices. Avron et 
al [|1^,0] showed that a generic Hamiltonian close to the degeneracy point can be expressed 



as 



Hap -^QijSl^^Si^p (67) 

where 5*^^ are the three 4x4 spin 3/2 matrices, and Qij is a 3 x 3 real symmetric traceless 
matrix. Since Qij has 5 real entries, and one needs to tune all of them to zero to get a 
degeneracy of two Kramers doublets, therefore the codimension of the degeneracy point is 5. 



In fact, this is nothing but a quadrupole spin Hamiltonian first used by Zee and other 
authors [0,0 to illustrate the concept of non-Abelian holonomy. 



This generic Hamiltonian for the symplectic class (p7| ) has a hidden 50(5) symmetry. 
Since Qij has exactly five entries, we can choose a orthonormal basis set Q^j, {a = 1,..5) 
satisfying Tr{Q°-Q^) = ^6°-'' and expand Qij as Qij = J2a^aQij- In this representation, the 
generic Hamiltonian ( p^ ) can be expressed as 

'Hafd = TT'oT'^i^ ; = Sl^^QijSi^p (68) 

At this point, the hidden 5*0(5) symmetry becomes manifest. Since the five F matrices obey 
a 50(5) Clifford algebra, this Hamiltonian is identical to the 50(5) spinor Hamiltonian (|TTp 



discussed in section III A 
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Therefore, we see that 5*0(5) symmetry arises naturally close to any generic degeneracy 
point of a time reversal invariant Hamiltonian with Kramers degeneracy. Both AF and dSC 
states in the high Tc problem are time reversal invariant, and their quasi-particle states 
form degenerate pairs. At the fermi liquid point, two of these pairs become additionally 
degenerate. Therefore 5*0(5) symmetry emerges naturally in this system. The non-Abelian 
holonomy of the SDW and BCS quasi-particles arises close to the fermi liquid degeneracy 
point, which is identified with a Yang monopole in the parameter space. 

III. EFFECTIVE ADIABATIC ACTION FOR THE SO (5) NON-LINEAR a 

MODEL 

A. Effective Action Arising from Berry's Phase 

Let us now consider the problem of a rigid 5*0(5) rotor interacting with a 5*0(5) spinor. 
We want to integrate out the spinor degrees of freedom and ask whether this procedure will 
generate any additional contributions to the action for the 5*0(5) rotator ( see Kuratsuji 



and lida and M. Stone |2J] for such discussion in the case of coupled 50(3) spinor and 
rotor ). 

Our system is described by the Hamiltonian 

n = no + ni (69) 

no = ^j:Ll (70) 

a<b 

Here Hq is the usual action of the 5*0(5) rigid rotor |l[ with Lab being the generators of 
rotation, and Hj is the Hamiltonian (|Tl|). The states of this composite system are defined 
as product states 

\^a[n\,na) = l^^aln]) ^\na) (71) 

where \na) are the states of the 5*0(5) vector and |\l'Q,[n]) are spinor states. We choose 
to be eigenstates of TCii^n, \1/), i.e. they diagonalize (|70D for a given orientation of n^. 
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=K|^«N) (72) 

The partition function is defined as 

ZiP)=J2{^4niO)],nM\e-^'^\\'^a[nm,nam (73) 

with summation going over all states of and \E'a- Let us perform time discretization and 
insert resolution of identity for the n-states. We make N intervals of length e with A^e = /3: 

Z{f3)= n Mnk) E (^ah(0)],n,(0)|e-^XiV-l))(n(iV-l)|...|^,[n(0)],n,(0)) (74) 

k=l •' n(0),a 

where dnijik) is an invariant measure of integration for the \na) states at interval k. When 
e is small we have 

{n{k)\e-^'\n{k - 1)) = {n{k)\e-^'''\n{k - i))e-WH"W)^ (75) 

and 

Z(/3) = E (^(0)|e-^«XA^ - l)){n{N - l)|e^^oXA^ - 2)) . . . (n(l)|e-^oXO)) x 

?l(0),Q 

X (^„[n(0)] le-^^^^^^^e-^^^"'^-^)^ . . . e-^^("^)^|^„[n(0)]) (76) 

The first line of the last equation is clearly the partition function of the 50(5) rotor. Fol- 
lowing Auerbach [^] we can write it as 

{n{k)\e-'^'>'\n{k - I)) = 

dfi{pk){n{kMk)){p{k)\e-'^'>'\n{k - 1)) 
= J dfi{pk) exp{ipa,k{na,k - na,k-i) - ^ ^(Pa'^fe ~ PbUaj'^e} 

a<b 

{na{k) - na{k - l)f 
= const X exp{-- — — — ^j-^ —e} (77) 

In the second line of equation ( fTBD we insert the resolution of identity at each r step 
i:/3j^ft.K])(^^Jn,]| = l,then 

Tna;na = [n(0)] | e-^("'^)^e-^("^-^)^ . . e-^^"^)^ | ^„ [^(0)] ) = 
n j:('^Mme-'^'^^^\^i>..AnN-i]){^b^.AnN-i]\ . . . \^a[nm) (78) 
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We now recall the adiabatic hypothesis and realize that in summation over intermediate 
spinor states in the last expression we only have to consider two states of the same energy 
as the initial state ( recall discussion of double degeneracy of spinor states in section [11 A| ). 
In the limit e ^ the last expression becomes the path ordered exponential 



Here 



Al,[n{t)] = -^{^i^4n]\^J^,[n]) 



(79) 



(80) 



and indices a, (3 take two values that correspond to the two degenerate states of the same 
energy as the initial state, but with the new orientation of the n-field ( new instantaneous 
eigenstates in the language of section |II A| ). Therefore we arrive at the same SU{2) gauge 
connection A"' as defined in equations (^6|) and (|1^). By combining equations (|76D and (|79D 
we see that this gauge field A"" defines a non-trivial non-Abelian topological term in the 
action for the SO (5) rotor 



Z{P) = jvurTr dr + i <f A''[n{T)]dna 



(81) 



In Chapter [ID we saw how the first term in this path integral may be written using two 



quaternions q = (qi, q2) 
1 / dua \ ^ 



(D.q)t(P,q) 



4\dT J 

Vrq= 9,q-qB(r) B(r) 



(82) 



As it turns out, the second term of equation (^Tj) can also be represented using quaternions 



iA°'[n{T)]dna qJ^S^q^fir = — B(r)(ir 
Therefore, the HP{1) representation of (|8T|) reads: 

Z{P) = JvBVciTr expl-J^^dr ((D,q)t(D,q) + B(r; 



(83) 



(84) 
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IV. PHYSICAL INTERPRETATION OF THE SU(2) HOLONOMY 



A. Holonomy on a Rung in the SO (5) Ladder Problem 

Having established the necessary mathematical framework, we are now in a position to 
discuss the physical application to the holonomy between SDW and BCS quasi-particles. 
As a warm-up exercise, let us first consider a two site problem ( see Figure §). 

— © — 



FIG. 3. Two site problem. Electrons on the upper site of the of the rung are created by cj. 
operators and dj^ creates electrons on the lower site. 

Recently, Scalapino, Hanke and one of us (SCZ) studied a 5*0(5) symmetric ladder 
model, and classified the SO (5) operator content for the two site problem. They introduced 
a 5*0 (5) spinor 



(c \ 
V4y 



(85) 



where c\ and (ij. operators create electrons on the upper and lower sites of the rung re- 
spectively. If we consider the 5*0(5) vector interaction on the same rung, and introduce a 
time-dependent Hubbard-Stratonovich field Uait) to decouple the interaction, we obtain the 
following effective fermion problem on a rung: 

n{t) = A n,{t)miTlp^> p (86) 

Our task is to find the ground state and low energy quasi-particle excitations for this problem 
and demonstrate the physical interpretation of the non-Abelian holonomy. 

This physical problem is identical to the mathematical problem we posed in section [11 C 



We can transform the problem into the "rotating frame" by decomposing according to 
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(p6|), with a Safi matrix satisfying (^) and (|30D. The resulting Hamiltonian ( pSD in the 
"rotating frame" is time independent and diagonal in the $o variable. $i and $3 have 
positive energy, while $2 and $4 have negative energy. According to the Dirac prescription, 
we fill the negative energy states and obtain the ground state \VL) in the rotating frame: 



(Notice that the definition of \VL) is different from reference [0). The corresponding elemen- 
tary excitations are given by 



with degenerate energy +A. If the vector points to n4 direction at t = 0, the ground state 
and four elementary excitations can be explicitly expressed in terms of electron operators 
Co- and dfj, and they are depicted in Figure ^. 



FIG. 4. Ground state and single particle excited states on a rung, when points in 714 direction. 

Within the adiabatic approximation, the Sap matrix is explicitly given by (pSf). The in- 
stantaneous ground state and elementary excitations in the original frame is simply obtained 
from ( |87D and (p8| ) by substituting $q, = 5'^^\l//3. Let us now imagine taking a path of na{t) 
on 5*^ which starts and returns to n^. Because of the non-Abelian holonomy discussed in 
section [11 C| , the Sap matrix only returns to itself up to a SU{2) transformation in the (1, 3) 
and (2, 4) subspaces. In particular, there exist cyclic paths on S'^ for which 



= = ; = $3|(^) = 



(87) 



(88) 
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Sit = 0) 



^ 1 ^ 



10 
10 
1 



Sit = T) 



'^0 10^ 



1 
10 
10 



(89) 



Under such a cyclic path, the quasi-particle state interchanges with ^I/slfi), while the 

quasi-particle state \E'2|f^) interchanges with \E'4|fi) ( see Figure |). 



adiabatic 
evolution 



adiabatic 

t evolution 

FIG. 5. Interchange of degenerate spinor states after cyclic adiabatic rotation of the superspin. 

Therefore, we obtained a simple physical interpretation of the S'f/(2) non-Abelian holon- 
omy in terms of interchange of certain degenerate pairs of quasi-particle states on a rung. It 
is important to notice that the degeneracy of these pairs is generic in any model with spin 
rotation and charge conjugation symmetry, and does not require any special fine tuning of 
parameters that is required to render the model S'0(5) symmetric in all sectors 0. For this 
reason, we believe that the non-Abelian holonomy we discussed in this work have general 
applicability beyond exactly 5*0(5) symmetric models. 

We notice that in the course of a cyclic adiabatic evolution, a SDW quasi-particle can 
be turned into a SDW quasi-hole (see Figure |^). In this sense, this phenomenon is very 
analogous to Andreev refiection at the metal/superconductor interfaces, where a metallic 
quasi-particle refiects back as a quasi-hole, and emits a Cooper pair into the superconductor. 
However, close inspection of the spin quantum numbers shows that a tt pair ll|,0], rather 
than a Cooper pair is emitted in our problem. Indeed, the states interchanged after cyclic 
adiabatic transport of the superspin differ from each other by addition/removal of a pair of 
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electrons in a triplet state ( a pair of spin up electrons for the first pair of degenerate states 
on figure and a pair of spin down electrons for the second pair of degenerate states on 
figure ^ ). Operators that create such triplet pairs are generators of the 50(5) symmetry 
on the rung (tt operators), they carry charge 2 and spin 1. By contrast, the operator 
that produces a d-wave Cooper pair creates two electrons on opposite sites in a singlet state, 
therefore it carries charge 2 and spin 0. 

The appearance of the vr operators as operators of the holonomy is not surprising and 
is related to the symmetry content of Hamiltonian (^6|). For a fixed direction of Ua the 
symmetry of this Hamiltonian is reduced from 5*0(5) to 50(4). For example when points 
toward n4, the 50(4) algebra that commutes with the Hamiltonian is formed by operators 
that rotate between rii, n2, n^, and but leave n4 invariant. Let Lab be the operators that 
rotate between and rih ( microscopically, such operators are given by Lab = ^ar^^\l/^, 
where F"* = i[F'^,F^], see p6|,|7| for details ), then the generators of the unbroken 50(4) 
symmetries are 

-ft'l = - {L21 - L35) -^2 = 1 (-^^13 - L25) K^ = - (L32 - L15) 

= - {L21 + L35) J'i — \ (-^13 + -^25) 2 ^^^'^ ~^ ^^^^ 

[Ki, Kj] = i eijkKk [Ji, Kj] = [Ji, Jj] = i tijkJk (90) 

We have chosen the generators in a form that makes transparent a property of 50(4) to 
be factorizable into SU{2) x SU{2). The eigenstates of our Hamiltonian have to form 
multiplets that transform as irreducible representations of this unbroken 50(4). But 50(4) 
multiplets also factorize into products of SU{2) multiplets, and by acting with operator 
= Tcl — in^ = — icjrfj or J+ = tt^ — iiTy = —icidi one can move between states related by 
holonomy. For example, states and ^al^^) form an 50(4) multiplet that is a doublet 

of the J algebra and a singlet of the K algebra, whereas the states $2!^) and $3!^) form 
an 50(4) multiplet that is a doublet of the K algebra and a singlet of the J algebra. Also 
since the two pairs belong to different multiplets holonomy does not mix them. 
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B. SU(2) Holonomy of BCS and SDW Quasi-particles 



Having illustrated the physical interpretation of the SU{2) holonomy in a 2 site problem, 
we are now ready to investigate the physical interpretation in terms of SDW and BCS 
quasi-particles. We start from the tight-binding Hamiltonian with S0{5) symmetric vector 
interaction. We write it in spinor representation following notations in 



pa 

+ E'Ap n^t) i^^^Tl^^^^df^ + ^UdaK^^rX^) (91) 

pa/3 



Here 



K = i^Ph Cpi, <Pn{p) C-p+Qt, Mp) C-p+qJ (92) 

is the 5*0(5) spinor with 4>n{p) = sgn{cospx — cospy) and prime over J2p' denotes summation 
over half of the Brillouin zone). When points in the 2, 3, 4 directions, this is nothing but 
the SDW mean field Hamiltonian. When Ua points in the 1,5 direction, it represents the 
BCS mean field Hamiltonian. We are interested here in the case where the na{t) field is time 
dependent. This situation would arise naturally when one performs a Hubbard Stratonovich 
decomposition of the interaction term, and na{t) is a dynamic mean field. As the na{t) field 
fluctuates in time, the quasi-particles "rotate" between SDW and BCS characters. Here 
we would like to consider the adiabatic limit of this fluctuation and identify the associated 
holonomy between SDW and BCS quasi-particles. 

Following the discussions in sections |1I Q , we decompose \E'pa as 

^pa(t) = Sa(iit)%f3it) (93) 

where the spinor rotation matrix Saj3 satisfies the same conditions (pQ]) and (|30D as before. 
Substituting this decomposition into (pl|), we find the Hamiltonian in the "rotating frame": 

^0 = E'ep(*!>a*P" - '^'L+Qa'f P+Qa) + E'^P {^U^%%+QP + ^Uq.^U^p) (94) 
pa pa/3 
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This Hamiltonian can be easily diagonalized by a 5*0(5) generalization of the Bogoliubov 
transformation: 



7^1 


= Up^pl + Vp%+Ql 


ft 

7-pT = 


Mp^'pS + ^^p'^'p+QS 




7pi 


= -Up^p2 + fp$p+Q2 


7-pi = 


-Up%4 + Wp$p+Q4 




7^1 


= ^^p'^'pl - Mp'^'p+Ql 


ct 

7-pT = 


^'p^pS - Mp'^'p+QS 




%i 


= Wp^'p2 + ^ip"^'p+Q2 


7-pi = 


yp^'p4 + Mp'^'p+Q4 


(95) 



where Vp = y |(1 - CpZ-Ep), Up = y|(1 + epZ-E'p) and Ep = ^t^ + . We obtain 

^0 = E (7pt7;. - llill.) (96) 

per 

and the ground state \VL) is obtained by filling the negative energy states: 

7^pjl^) = = (97) 

Excited states 

7±U|f^) 7±p.l^^) (98) 

are degenerate with energy Ep. If at t = the superspin points towards n^, these 
quasiparticle states are readily expressed in terms of electron creation/annihilation operators 

7±pal^) = (^pC±p<x - MpC±p+Q^) \n) (99) 

Elementary excitations in the original frame are constructed as $pQ, = S'^^(t) \E'p/3 where in 
the adiabatic approximation matrix Sais is taken from (pSf ). 



Now we imagine that na{t) makes a closed adiabatic rotation as discussed in section |IV A . 
It starts from and returns at t = T to the same direction. However, the matrix 5* at time 
T does not have to return to unity and it may come back in the form given in (|89D. Equation 



(p3|) then instantly tells us that after such cyclic rotation of the superspin the quasiparticle 
state 7p||f2) interchanges with 'j^^-^\Q), and 7p||fi) interchanges with j^^^lQ). This means 
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that after superspin completes a revolution, quasiparticle in the conductance band returns 
as a missing quasiparticle in the valence band with the same spin but opposite momentum. 
SDW quasiparticle and SDW quasihole have been interchanged! What has been emitted in 
such process is an object with spin 1 and charge 2 



7'p<x7L = Mp) C-p+QaCpa + ■■■ (100) 



When summed over all p's the terms in . . . disappear and we get one of generators of 5*0(5) 
algebra, the vr operator . So we obtained again that adiabatic transport of the superspin 
leads to a S0{5) generalization of Andreev reflection, in which a vr particle has been emitted. 



Interpretation of this result using unbroken SO (4) symmetry of Hamiltonian (^) is 
straightforward. When points along n^, the important generators of this unbroken sym- 
metry ( see discussion after equation (pO| ) ) are 

= i(L2i - L35 + iLi3 - iL25) 

p p 

J+ = ^(-^21 + L35 + zLi3 + iLas) 

= E<^-(p)c-p+QTC-pT = E 7-pt7pT (101) 
p p 

These are n operators that move us between pairs of states that form multiplets of the 

unbroken SO {4). 

7pil^)-7-Ul^) 
7-pJfi) - 7?J^^) 

7-pT|fi)^7Sl^) (102) 

Arguments presented in this Section clarify the origin of the non-Abelian holonomy of 
spinors in the 5*0(5) theory. It comes from a non-zero overlap between SDW and BCS 
quasiparticle states. 
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V. FERMIONS AND THE SO(5) NON-LINEAR a MODEL 



We can now put our knowledge to a good use by considering how fermions can be added 
to the 5*0(5) Non-Linear a ModeL We start with a microscopic Hamihonian with SO{b) 



symmetric pseudo- vector interaction ( spinor notations of [36] are used again ): 



n = -tY, ^>\x)'^{x + 5,) + At (x) A J 



X] 



^a{x) = E(-)'^^(a^)^r'^*(a; + 6,) (103) 

where is a translation by one lattice constant in direction fi, and (— )^ = 1 when fi = x 
and —1 when fi = y. In momentum space such interaction corresponds to 

H.„i = \/5:At(q)A,(q) 
q 

A„(q) = 2 Wk^k+qi?r«^_k (104) 

k 

with = cos{kx) — cos{ky). The mean field of this Hamiltonian in direction Ai readily 
reproduces the usual d-wave superconductivity. Antiferromagnetic phase is less conventional 
since the S0{5) symmetry dictates it to have nodes in the energy gap, with no sign change 
of the order parameter (see for a more detailed discussion). 



We can use Hubbard-Stratonovich transformation to write Hamiltonian (|103D as 



n = -tY: v&nx)v]/(x + 6,) + At(x)A„(x) 

+ J2 ^li^) + '^m) + h.c. (105) 

A simple Hartree-Fock saddle point for this Hamiltonian has a disadvantage that it does not 
reproduce the non-linear a model directly. But we can do better than the mean-field if we 
explicitly separate degrees of freedom that contain Goldstone bosons of the SO (5) rotations 
(see [0,0 for analysis of the 5*0(3) symmetry in the Hubbard model and [^ for analysis 
of U{1) symmetry and nodal quasiparticles in d-wave superconductors ). It is important to 
note that in the long- wavelength hmit the fields Aa{x) are real ( At(q = 0) = Aa(q = 0) 
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), so we can use matrix S from Section |II C| to rotate all spinors locally to point in one 
direction, i.e. we introduce new spinors as 



where rotation matrix S{x) satisfies two conditions 

Aa{x)S\x)T''S{x) = r^A(x) 
S\x)S{x) = 1 

with A^(x) = Y,a'^i{^)^a{x). One can easily prove for S the identity 

S^RS = R 

and use it to rewrite the first equation of ( |107| ) as 

Aa{x)S^{x)Rr''S{x) = Rr^A{x) 



(106) 



(107) 



(108) 



(109) 



Neglecting the high energy amplitude fluctuations of A{x) we can express Hamiltonian (|105|) 



as 



Umf = -tYl + 5^) + AY,{-Yi>{x)RT^^{x + 5^) + h.c. 



x,5u 



Hi = -tJ2 ^^x) (s^{x)S{x + 5^) - l) ^{x + 6^) + h.c 



^2 = {S^{x)RT''S{x + 6^)Aa{x) - RT^a) ^{x + 6^) + h.c. 

It is helpful to define \^(a;) = G^{x). Then one can use 

GRV^G = 

to write the mean-field Hamiltonian as 



(110) 



^ i(7y ^ 



V 



—io,, 



(111) 



71+ = Y^ {ek($kl^'kl + ^'k2'^'k2) + AWk(<lfei<l_k2 - <i>k2$-kl + '^L^'-kl - '^'kl'^'-k2)} 
k 

^- = E {ek(<i'L*k3 + ^L*k4) - A^^;k($k3$-k4 - $k4$-k3 + <^L^tk3 - ^L^-k4)} (112) 
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with ek = -2t{coskx + cosky). We introduce = {l>ki, $k2, -$Lk2' "^-ki} ^-^^ 
{$k3, ^k4, -^'-k4' "^-ks) w"te equation (|112|) as 



(113) 



where we used r matrices 



T"3 



^ ^ 

-i 



(114) 



and later we will also use a\ matrices defined as 



cyx 



(115) 



We expand (|113|) around Ki = (|, |) and use cki+q ~ —2t{qx + qy), WKi+q ^ Qy — Qx 
with analogous expressions for K2 = (— f , f )• If we now define $±^K,+q = 0±iq "we can write 
the mean-field Hamiltonian as 



q,r=± 

+ H 0l2q(2t(-gx + gy)T-3 + rA(g^. + gy)ri)0,.2q 



(116) 



q,r=± 



It is convenient to introduce new coordinates xi_2 = and velocities vi = 2y/2t, V2 
2\/2A, then from ( |116|) we obtain free Dirac Hamiltonian in 2 + 1 dimension 

d \ , 



+ rv2{-tTi, 



_d_ 

dx-} 



_d_ 

dxi 



} 



(117) 



Let us now look at Tii. From equation (0) of Section [11 C| we know that in the subspace 
of components 1 and 2 we have S^d^^S = iA^^^crx, and in the subspace of components 3 
and 4: S'^df.S = iA^^^^ax- Here A^^^ = ELi ^^±^^^nb with Af given by equation (|19|) (see 
also (0)). Other elements of S'^d^S are zero in the adiabatic approximation. Therefore, Tii 
can be expressed as 
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ni = -tJ2 ^\x)G [s\x)S{x + 6f,) - l) G<l(x + 5p 



-tE<|t(a:) 



\ 



lAl^^ax 
lA^L 



\ 



(118) 



Taking the $+,k components of the last equation first, we have 



kk'q/i 



= 2tY: '$;,k^t,/.(q)^Ar3$+,k'Sm(A;;) 5(k - k' + q) 

kk'q/i 

Expanding the last expression around Ki and K2 we obtain 



iAlJq)ax 



v 



$k'2 



sin{k')6{\i - k' + q) 



(119) 



ni+ = 2t (^+,..(q)+^+,,(q)XKi+q,^Ar3$+,K,+q,5(qi-q2 + q) 
qiq2q 

+ 2t E (-^+,x(q)+^+,,(q)XK.+q,^AT-3$+,K.+q.5(qi-q2 + q) 
qiq2q 

= vi j dx [4>\^{x){-iT:i){iAX^^{x)ax)<p+i{x) + 4>\^{x){-iT^){tAX^^{x)dx)<p+2{x)] (120) 



where ^ = {A\^^ + y)/v/2 and 2 = {^X,x - ^X,y)/^- Analogously, for the $_,k 
components of Tii we have 

ni- = vi j dx{0li(x)(-zr3)(z^^,i(x)aA)0-i(a:) + 0^_2(a;)(-^r3)(U^,2(a;)^A)0-2(x)} (121) 



Similar manipulations may be performed for 7-^2- One uses equations (|107|) , ( p.09| ), and 
( |1 1 1|) to prove the identity 

G ( S'^{x)RT''S{x + 5^) /V{x) - i?r^A ) G = A GRT^G G ( 5^(x)^(x + 5^)-l)G 



( 



A 



V 



A^ 



(122) 



and then follows the steps that lead to ( |12(J| ) and ( |121D . After straightforward manipulations 
we obtain 
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kk'q/x 



Combining all the pieces together we obtain 

d 



dxidx2{ vi 4>li{x){-in)[-^ + iA^.^iax](l)r 



Ax 



d 
d 
d 

+ rv2 0l.2(a;)(-«ri)[— + iA^.^ax]<pr2{x) } (124) 

OXi 

It is interesting to notice that we have a two-dimensional representation of the Dirac matrices 
for our Dirac fermions. So each of the ipri and 0^2 are chiral fermions, but parity is not broken 



since parity transformation (t—^t,y-^y,x—^ —x pO|,^ ) simply interchanges 0^1 and 



<Pr2- 

Combining equation ( |124| ) with the HP{1) representation of the 5*0(5) non-linear a 
model (|55D, we see that the SO (5) theory including both fermions and bosonic superspin 
variables can be completely formulated as a 2 + 1 dimensional relativistic SU (2) gauge field 
theory, with two bosonic "Higgs fields" qi and q2 and four flavors of Dirac fermions (pra 
(r = ± and a = 1, 2) [|^. The gauge field does not have a kinetic term J-"^^ in the ordered 
phase. However such a term would be generated in the quantum or thermal disordered phase 
]. This complete formulation of the SO {5) theory is a central result of this work. 



VI. CONCLUSIONS 

Our work showed that the "square root" of the 5*0(5) theory, i.e. the spinor sector, 
has a fascinatingly rich internal structure. From the point of view of the 50(5) theory, the 
fermionic quasi-particles of the high Tc superconductors are still the ordinary SDW and BCS 
quasi-particles. The novel aspect of this system lies in the interplay between them. Our 
paper is a tale about the intimate relationship between these two quasi-particles. The path 
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leading towards their union travels through some of the most beautiful areas of modern 
mathematics, with 5*0(5) symmetry being a unifying central theme. This mathematical 
structure enables us to formulate the complete 5*0(5) theory, including both fermionic and 
superspin variables. 

In the following, we shall discuss some possible extensions of our work. 

5*0(5) symmetry breaking: In order for these mathematical ideas to apply to the real 
systems, the effects of explicit 5*0(5) symmetry breaking has to be carefully addressed. 
As noticed in reference [^, the degeneracy of the 5*0(5) spinor multiplet is much more 



robust than the vector multiplet, therefore, the results obtained in the work could have 
more general applicability beyond exact 5*0(5) symmetric models. In a model with 5*0(5) 
symmetry breaking, the superspin is no longer confined to move on 5*^, but will trace out 
a general trajectory in the five dimensional space. However, as long as it does not pass 
the degeneracy point at the origin, the topological nature of the holonomy ensures that the 
results would still apply in this case. 

Non-Abelian Bohm Aharonov effect: In principle, one could construct AF and SC het- 
erostructures and manipulate the direction of the superspin by magnetic fields and currents. 
Since the 5*f/(2) gauge field is explicitly determined by the superspin direction, one could 
construct regions with finite SU (2) flux, split the quasi-particle beams around it and observe 
their interference pattern. The Abelian Bohm-Aharonov interference produces modulation 
of the particle density, while the non- Abelian Bohm-Aharonov interference produces modu- 
lation of the mixing ratio of the particles belonging to the 5*f/(2) doublet. 

Non-trivial fermion numbers of topological defects: Recently, many authors investigated 
various types of topological defects in the 5*0(5) theory |1T0|,^, |4^JTT| . These defects usually 
involve non-trivial spatial variations of the superspin direction, their associated SU (2) gauge 



field could give rise to non-trivial fermion numbers around the defects . 

50(5) Andreev reflection: In Section |V|we have seen that propagation of quasiparticles 
through regions with nonuniform direction of the superspin may lead to 50(5) Andreev 
reflection, when SDW quasiparticles in the conduction band turn into SDW quasiholes in 
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the valence band or vice versa. This process is similar to Andreev reflection at the super- 
conductor/normal metal interfaces, with an important difference that the emitted particle is 
not a Cooper pair, but a n pair, a triplet two particle excitation. Such processes could lead 
to novel phenomena in superconducting/ antiferromagnetic heterostructures: the appearance 
of resonant tunneling as in or the possibility of new bound states in heterostructures 
and around topological defects of 5*0(5). Our work gives a natural 5*0(5) generalization of 
the Bogoliubov deGennes formalism to treat general bound states in all these cases, with 



Sommerfeld quantization condition for the existence of the bound state becoming a 2 x 2 
matrix equation due to an SU{2) holonomy of the 50(5) spinors. 

Single particle spectra in the pseudogap regime: Within the 5*0(5) theory, the pseudogap 
regime is interpreted as the fluctuation regime of the 5*0(5) superspin vector. The quasi- 
particles in this regime have fluctuating SDW/BCS characters, and can be naturally treated 
within the flnite temperature formalism of Dirac fermions coupled to fluctuating 5*f/(2) 
gauge fields. Connections to the photoemission experiments in the pseudogap regime could 
be made. 

Relationship to other works: The 5*f/(2) holonomy of BCS and SDW quasi-particles dis- 
cussed in this paper bears some formal resemblance to the Affleck-Marston SU (2) symmetry 



48| in the flux phase and the SU (2) gauge theory recently formulated by Lee, Nagaosa, Ng 



and Wen However, these two approaches are obviously motivated by very different 



physical pictures and have entirely different physical meanings. Nevertheless, it would be 
useful to explore their formal connections. 

Recently, Balents, Fisher and Nayak studied the evolution of the gapless fermion 



spectrum around the d wave nodes from the SC state to the insulating state. If the insulating 
state in question is a AF state, the formalism developed in this work might be related to 
some of their considerations. 

High Tc superconductivity involve strong electron correlations which usually makes per- 
turbative calculations difficult to carry out. On the other hand, topological effects are robust 
and have general validity even in the strong interaction regime. Our work attempts to lay 
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down a mathematical foundation for future studies of topological effects within the 5*0(5) 
theory. We hope that Nature made use of these elegant mathematical concepts in the high 
Tc superconductors. 

We would like to thank Prof. A. Auerbach, E. Fradkin and A. Zee for useful discussions. 
This work is supported by the NSF under grant numbers DMR-9400372 and DMR-9522915. 
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